function [coords] = quantify(coords, deltaR, deltaPhi)
[phi rho] = cart2pol(coords(:,1), coords(:,2));
rho = deltaR*floor(rho ./ deltaR);
phi = deltaPhi*floor(phi ./ deltaPhi);
[coords(:,1) coords(:,2)] = pol2cart(phi, rho);
end